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I. INTRODUCTION 

Squeezed states of light offer possibilities of improving the performance of optical devices 
since they can reduce fluctuations in one quadrature below the level associated with the 
vacuum states. Such a situation is relevant for the optical communication networks as 
well as for many optical devices. For instance, recently such light has been used in a 

n n 

power recycled interferometer [1| and in a phase-modulated signal-recycled interferometer [2|] 
aiming to improving significantly the sensitivity of these devices. It has been shown that this 
light can be used to tune the resonant frequency of the cavity without actually moving the 
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signal recycling mirror or changing the bandwidth of the interferometer without substantially 
decreasing the sensitivity at the resonant frequency J2|. Also squeezed light has been applied 
in quantum information theory, e.g. in a quantum teleportation 3], cryptography 4J, [5| and 
dense coding [6|]. In this respect the security in a quantum cryptography 5| relies on the 
uncertainty relation for field quadrature components of these states. Further, experiments 
on quantum teleportation have been successfully performed by means of two-mode squeezed 
vacuum states [7|. A variety of methods have been proposed for the generation of squeezed 
states of the electromagnetic field and several of them have been realized experimentally 
[Si, |9] . For more complete information about squeezed light the reader can consult review 
papers lOj. 

As is well known the ideal single-mode squeezed light (with possible 100% squeezing) 
may be generated using a degenerate parametric amplifier (DPA) with classical pump. This 
is the most effective system in nonlinear optics in which such strong squeezing was obtained. 
DPA is the two-photon process where the signal and idler modes are identical in this case. 
Such a property manifests [111] itself in the large scale oscillations in the photon-number 
distribution [l2[ and in the bifurcation in phase distribution [13]. In other words, the 
occurrence of squeezing in the quadrature variances does not need to be accompanied by 
oscillations in the photon-number distribution or bifurcation in the phase distribution. The 
best example are the binomial states |l4| . On the other hand, the nondegenerate parametric 
amplifier (NDPA) can produce ideal squeezing in the compound modes even if the single 
mode squeezing does not appear 15j |. 

In this article we suggest a compound system-nonlinear three-mode asymmetric direc- 
tional coupler which is able to provide both strong single mode squeezing and compound two 
mode squeezing under certain conditions and for specific modes. This means that the present 
system can operate simultaneously as degenerate and nondegenerate amplifiers. Moreover, 
the single mode squeezing can be confirmed in the behaviour of both the photon-number 
distribution and phase distribution. In general, the nonlinear directional coupler consists 
of two or more parallel optical waveguides fabricated from nonlinear material. Both waveg- 
uides are placed close enough to permit flux-dependent transfer of energy between them as 
a consequence of evanescent waves involved in the interaction [16j. This flux transfer can 
be controlled by the device design and the input flux. Further, the nonlinear directional 
coupler has been experimentally implemented , e.g. using semiconductor/glass composite 



17|, erbium doped fibers 18|, ll9| and a crystalline organic semiconductor 20] . Furthermore, 



the directional couplers are of the technological interest since they are currently a topic of 
major interest in integrated optics for ultra-high-speed optical signal processing and optical 



switching 



21] as well as in the generation and transmission of nonclassical light. In the view 



of these reasons it is certainly useful to make concrete efforts to study the actual possibility 
to generate and to transfer squeezed light in this system. 

It is important to mention various works investigating similar systems, e.g nonlinear 
asymmetric directional couplers which consist of a linear waveguide and a nonlinear waveg- 
uide operating by second harmonic generation 22H24J or of a nonlinear waveguide with the 
core of Kerr-like medium 25[ ; it has been shown that such systems can be used to reduce the 
necessary switching power (for the switching of a signal controlled by a strong pump) 26] 
as well as they are useful for constructing a band-pass power-filters or a band-reject power 
filte r |25] . For a review of the role of quantum statistical properties in nonlinear couplers, 
see [22 1 • 

In the present work, we show that a nonlinear asymmetric directional coupler composed 
of a linear waveguide and a nonlinear waveguide operating by nondegenerate parametric 
amplification is an effective source of squeezed light. This will be done as follows: In section 
2 we describe the model under discussion and give the solution of the equations of motion. 
In section 3 we demonstrate the single-mode quasidistribution functions. Sections 4 and 5 
are devoted to the photon-number distribution and phase distribution, respectively. Finally 
we summarize the main conclusions in section 6. 



II. EQUATIONS OF MOTION 



We consider a device (nonlinear asymmetric directional coupler) as outline in Fig. 1. It 

can be described by the Hamiltonian 
H 3 



T = C w i 



h 



oJ-Oj] 



iAi {aia2 exp[i(u;i + uj-zjt] — h.c.} 



i=i 



-MA 2 I aia\exp[i(ui — u 3 )t] — h.c. > + i\ s < c^a^exp^u^ — u> 3 )t] — h.c. > , (1) 

where %(aj), j = 1, 2, 3 are the annihilation (creation) operators designated to the signal, 
idler and linear modes, respectively, Ai and Xj, j = 2,3 are the corresponding nonlinear 
and linear coupling constants; u)j are the natural frequencies of oscillations of the uncoupled 



modes and h.c. is the Hermitian conjugate. The linear exchange between the two waveguides 
establishes through evanescent wave provided that u\ ~ u 2 — 0J3. Moreover, we treat the 
problem of propagation in the Hamiltonian formalism neglecting dispersion. Thus if all 
waves are propagating with the same velocity, time t and space z relate by the velocity of 
propagation v, z — vt. For further details of quantum description of propagation we refer 
(and references therein). On the other hand, as we assume one-passage 



the reader to 



propagation, losses in the beams can be neglected. In general they can be described in 
the standard quantum way in the form of interaction of light beams with reservoirs, as for 
instance described in 29] . However, one can note that nonclassical properties of light beams 
degrade by the effect of the reservoir, less by the damping which may be by several orders 
weaker than the nonlinear coupling, more by the influence of non-zero mean numbers of 
reservoir oscillators. 
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Fig.l: Scheme of realization of interaction in (1) using a nonlinear asymmetric directional 
coupler which is composed of two optical waveguides fabricated from first-order (x^ 1 ' ) and 
second-order (x ) materials, where x designate susceptibility. Signal mode 1 and idler mode 
2 propagate in the first waveguide and linear mode 3 in the second waveguide. The interaction 
between the signal and idler modes is established by strong pump coherent light, which is 
not indicated in figure, with the coupling constant \\. The interactions between the linear 
mode and the signal and idler are established linearly with the coupling constants X 2 and 
A3, respectively. The beams are described by the photon annihilation operators as indicated; 
z = vt is the interaction length and we assume that all beams have the same velocity v 
and the length of the waveguides is L. Outgoing fields are examined as single or compound 



modes by means of homodyne detection to observe squeezing of vacuum fluctuations, or by 
means of a set of photodetectors to measure photon correlations, photon antibunching and 
sub-Poissonian photon statistics in the standard ways. 

The dynamics of the system is described by the Heisenberg equations of motion which, 
using the slowly varying forms (Sj = Aj exp(— iuijt), j = 1,2,3) to eliminate free oscilla- 
tions so that only the slow dynamics due to the coupling between modes appear explicitly, 



read ^ 
gL4i 

~dT 

dA 2 

dt 

dA 3 

dt 



-Aii 2 - A 2 i 3 , 
-Aiil - A 3 i 3 , 



x 2 A 1 + \ 3 A 2 . (2) 

These are three equations with their Hermitian conjugates forming a closed system which 
can be solved easily by the Laplace transformation, restricting ourselves to the case A 2 = 
A3 = ■%, to avoid the complexity in the calculations; this means that we consider stronger 
nonlinearity in the first waveguide compared with the linear exchange between waveguides, 
thereby having the solution 

Mt) = ai(0)A(t) + ol(0)/ 2 (t) - a 2 (o)/ 3 (t) 

-al(0)f 4 (t) - a 3 (0)/ 5 (t) - a 3 (0)/ 6 (t), 



M*) = 02(0)01 (*) + al(0)g 2 (t) - oi(0)^ 3 (t) 
-o t 1 (0)^ 4 (*) - a 3 (0)g 5 (t) - a 3 (0)(? 6 (t), 



Mt) = Os(0)M<) + 4(0)/i 2 (t) + a 2 (0)h 3 (t) 

+a\(0)h4t) + Oi(0)/i 5 (t) + d\(0)h 6 (t), (3) 

where the time-dependent coefficients, i.e. fj{t),gj(t), hj(t), including all information about 



the system are 

hAt) = \ 



Mt) = 



cosh(Ait) ± cosh( — ) cos(kt) ± —= sinh( — ) sin(H) 
sinh(Ait) =F sinh( — t) cos(kt) =F —= cosh( — ) sin(kt) 



2 At 

-cosh(^-)sin(H), f e (t) 



-sinh(^-)sin(H), 



(4) 



h±(t) = cosh( ) cos(kt) — sinh( ) sin(kt), 



h 2 {t) = — sinh( — ) cos(kt) -\ — -= cosh( — ) sin(kt), 
2 y3 2 

-cosh(— )sin(fct), h 4 , 6 {t) = -J - sinh(— ) sin(fct). (5) 

where k = ^Ai and the expressions for gj(t) are the same as those of fj(t). 

In fact, the nature of the solution can show how the coupler does work. To be more spe- 
cific, the time- dependent coefficients contain both trigonometric and hyperbolic functions. 
Consequently, the propagating beams inside the coupler can be amplified as well as switched 
between waveguides, i.e. between modes, in the course of time. 

On the basis of the well known commutation rules for boson operators, the following 
relations can be proved for the time-dependent coefficients 

fi(t) - f 2 (t) + m - fiif) + m - m = 1, 

fi(t)gi(t) - f2(t)g 3 (t) + h(t)g 2 {t) - U(t) 9l {t) - f 5 {t)g 6 (t) + f 6 {t)g 5 (t) = 0, 
fi(t)g 3 {t) - f2(t)g,(t) + f 3 (t)gi(t) - U(t)g 2 {t) - h{t)g 5 (t) + f G {t)g 6 (t) = 0. (6) 

The remaining relations can be obtained from (6) by means of the following transformations 

(fi(t)j 2 (t),Mt)j 4 (t),h(t),f 6 (t)) 

< — y (-#»(*), -04(*),-</i(*), -02(f), <fe(t),06(*)J 

^ (h 5 (t),h 6 (t), -h 3 (t), -h A (t), -h^t), -/i 2 (£)) . (7) 

Based on the results of the present section, we can study the quantum properties of the 

evolution of different modes in the model when they are initially prepared in coherent states. 

This will be demonstrated in the following sections. 

III. QUASIPROBABILITY FUNCTIONS 

Evaluation of various time-dependent mode observable is most conveniently achieved 
with the aid of corresponding time- dependent characteristic functions, their normal, anti- 
normal and symmetric forms, and the Fourier transforms of these characteristic functions 
(quasiprobability functions). All of these are related to the density matrix which provides 
a complete statistical description of the system. There are three types of quasiprobability 
functions: Wigner W-, Glauber P-, and Husimi Q-functions. These functions could be used 
also as crucial to describe the nonclassical effects of the system, e.g. one can employ the 
negative values of jy-function, stretching of Q-function and high singularities in P-function. 
Furthermore, these functions are now accessible from measurements [301 . 



7 

In the following we consider phase space distributions for the single-mode when all modes 
are initially prepared in coherent states before entering the coupler. 

We start from knowledge of the single-mode s-parameterized characteristic function which 
is suitable to describe the quantum statistics of light and for the jth mode it is defined by 

C,-(C, t, s) = Tr {p(0) exp [(A](t) - CMt) + |lC| 2 ] } , (8) 

where p(0) is the initial density matrix for the system under consideration, j takes on the val- 
ues 1, 2, 3 corresponding to considered mode and s takes on values 1, and —1 corresponding 
to normally, symmetrically and antinormally ordered characteristic functions, respectively. 
Assume that the modes are initially uncorrelated and are described by the density operator 

3 

p(o)=ni«)jj(«i' ( 9 ) 

where \a)j are coherent states. So that equations (8) and (9) lead to the single- mode s- 
parameterized characteristic function which for the signal mode takes the form 

CMC, s,t) = exp ji|C| 2 [s - \m(t)\] + \V2(t)\(( 2 + C 2 + [al(t)C - «i(*)C*]} , (io) 

where a\(t) is the mean value of the operator A\{t) in the coherent state and rjj(t), j = 1,2, 
are given by 

m(t) = m + m + m + m + m + m, 

r*(t) = fi(t)f 2 (t) + h{t)h(t) + h{t)U{t). (11) 

The single-mode s-parameterized quasiprobablity function is defined as 

Wj(P,t,s) = ^ J C j ((,t,s)exp(PC - /3*C)rf 2 C, (12) 

where Cj((,t,s) is the s-parameterized single-mode characteristic function and when s = 
1,0, —1 relation (12) gives P-, W- and Q-functions, respectively. 

Now the s-parameterized single-mode quasiprobability functions of the signal mode, 
which can be derived from (10) and (12), are 

WiOW) = — 2 

iry/iA+tf) - s][A-(t) - s] 

nm-mi m+mn (l3) 

XeXP \2[A.(t)-s] ~ 2[A + (t)-s]j> (13) 

where l±(t) = a±(t) — (3, a\(t) has the same meaning as before and 

A ± (t) = [h(t) ± f 2 {t)f + \h{t) ± U{t)f + Mt) ± /e(t)] 2 ; (14) 

the other expressions for the idler and linear modes can be obtained from (13) and (14) 

by making use of the transformations (7); we except here the case s = 1, which will be 

discussed at the end of this section. Now let us start our investigation by demonstrating 

the behaviour of the W^-function, i.e. setting s = in (13). In the language of mechanical 
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analogy of a harmonic oscillator with dynamical conjugate variables x and y, i.e /3 — x + iy, 
the VF-function can be written in the form 

where 

((Al(t)) 2 ) = ^, ((Al>(t)) 2 ) = ^1, (16) 

where a\{t) = a x (t) + ia y (t); ((AX(t)) 2 ) and ((A7(t)) 2 ) are x- and y-quadrature vari- 
ances, respectively. As known these quadratures relate to the conjugate electric and mag- 
netic field operators of the electromagnetic field and may be used as a measure of squeez- 
ing phenomenon. More illustratively, the two quadrature operators of the jth single- 
mode are defined through X(t) = \[Aj(t) + A](t)\, Y(t) = j-[Aj(t) - A){t)}, where 
[X(t),Y(t)] = | and then the uncertainty-relation reads ((AX(t)) 2 )((AY(t)) 2 ) > i where, 
e.g. {(AX (t)) 2 ) = {(X(t)) 2 ) — (X(t)) 2 . Therefore, we can say that the jth mode is squeezed 
if 4((AI(t)) 2 ) - 1 < or 4((Ay(t)) 2 ) - 1 < 0. Now from (14) and (16) together with (5) one 
can easily prove that the quadrature variances of the linear mode (single-mode squeezing) 
are 

((AX(t)) 2 ) = 1 |^sin 2 (H) + [cos(H) + -L s in(H)] 2 } exp(-A 1 t), (17) 

((A7(i)) 2 ) = i |^sin 2 (jfet) + [cos(H) - ^sin(H)] 2 |exp(Ait), (18) 

where k has the same meaning as before. Thus squeezing can be achieved in the linear 
waveguide in the X-quadrature. Moreover, it is clear that squeezing values become more 
pronounced for a large interaction time, i.e. for long length of the coupler L. This behaviour 
shows that changing the power of the linear interaction it is possible to transfer nonclassical 
properties, which are generated in the nonlinear waveguide, to the linear signal mode. For 

completeness, the uncertainty relation reads 

1 16 

((AX(t)) 2 )((AY(t)) 2 ) = -[l + ^sin 4 (H)]. (19) 

This formula reveals that the minimum- uncertainty relation holds only when t — t s — ^S^- 
where m is a positive integer and we call t s a squeezed time. In this case the device provides 
the squeezed coherent light in the linear mode with squeeze parameter r = 2 ^ L . It is 
important to mention that at t — t s the two waveguides become completely independent (cf. 
equations (3)- (5)) and the signal and idler modes can display perfect two-mode squeezing, 
which one can easily check with the help of two-mode quadrature variances. This can be 
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explained as follows: A portion of the energy always remains within the guide into which 
the field was initially injected. This energy grows in the course of time until t = t s , when 
the two guides are completely independent and hence the modes are trapped in their own 
guides. So that the self-interaction of the linear mode can attribute the well known single- 
mode squeezed light in the linear waveguide having its origin in the nonlinear waveguide 
in which the nonlinearity produces two- mode squeezing. For later times t ^ t s , the device 
generates single-mode squeezed light only in the linear waveguide where the corresponding 
VF-function may be broader than that of the squeezed coherent states. This situation is close 
to that of a two-photon absorber (a two-photon absorption by a reservoir of two-level atoms 
from a single mode of the electromagnetic field |31j) where a squeezed state, which is not 
a minimum-uncertainty state, has been generated [32| . The origin of squeezing of initially 
unsqueezed light interacting with two-photon absorbers is that the squeezing is generated 



by simple quantum superposition of states of light [31 1 . 

We proceed in our discussion by focusing our attention on the P-function. It is known 
that the correspondence between quantum and classical theories can be established with 
the use of this function. However, P-representation does not have all the properties of 
a classical distribution function, especially for quantum fields. To be more specific, light 
fields for which the P-representation is not well-behaved distribution (in most processes 
involving an interaction at least for some values of interaction time, including the process 
under consideration) can exhibit nonclassical features. We have shown earlier that this 
model is able to provide squeezed light and this should be reflected in the behaviour of the 
P-representation, i.e. setting s = 1 in (13). The significant example for this situation is the 
behaviour of the linear mode. For this case, we can show that the single-mode quadrature 
squeezing is established provided that A + (t) — 1 < and A_(£) — 1 > 0. It is evident that 
the P-function is not well-behaved function in this case and this is the indication of the 
nonclassical fields. Indeed, if we look at the model under consideration as a competition of 
parametric processes, we can find that this result is in marked contrast to the most results 



occurred in the literature for three interacting modes |33l. I34J . where the delta function of the 
initial P-representation of coherent light for a single mode becomes well-behaved distribution 
during the interaction. 

In the following sections we use the phase space distribution functions to study the 
photon-number distribution and phase distribution for the system under discussion. 
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IV. PHOTON-NUMBER DISTRIBUTION 



Photon-number distribution, i.e. the probability of finding n\ photons in the signal mode 
at time t ; can be obtained in the photodetection process, and is determined by means of the 
relation [33 1 



P(m,t) = (ni|pi(t)|m) 



\P 



2m 



ni } . 



-W(P,s = l,t)exp( 



Vfl 



(20) 



substituting (13) into (20) and carrying out the integration, we get the photon-number 

ai(t) - a\{t)} 2 [ai(t) + a${t)] 2 



distribution of the signal mode as 
P{ni,t) 



y/lA+it) + l][A-(t) 



+ 1] 

- n±—r 



: CXp 



2[A_(t) + l] 2[A+(t) + l] 




ai{t)-a\{t)f 
Al(t)-1 
where L*(x) is the associated Laguerre polynomials defined as 

. kt > A T(n + k + l)(-x) m 



L*Jx) 



m=0 
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m)\T(m + k + l)m\ 



(21a) 
(216) 



while T is the Gamma function. 



It is interesting to compare this distribution with the corresponding Poisson distribution 
P(ni,t) = Mll^ex P (-(fi 1 (t)}), (22) 

Till 

which corresponds to fully coherent field with the same mean photon number (hi(t)). The 
expressions for the idler and linear modes should be obtained from (21a) using the transfor- 
mations (7). 

In Fig. 2, we have plotted P(ni,t = 1.5) against n\ for the signal mode (solid curve) for 
ctj = 3exp(i|) and Ai = 0.3. For the sake of comparison, the corresponding photon dis- 
tribution for coherent field, (22), is shown by dashed curve. For all numerical calculations 
Yum=oP( n i>t) = 1 with good accuracy. We noted for the signal mode that the photon- 
number distribution exhibits always one-peak structure. From this figure we see that the 
behaviour is rather super-Poissonian as a result of quantum fluctuations because the solid 
curve is always broader than the corresponding dashed curve. So that, in general the Pois- 
sonian light evolves in the nonlinear waveguide as super-Poissonian light. On the other 
hand, as we have realized earlier from the behaviour of VF-function that squeezed states 
can be generated. However, these states are exhibiting oscillating photon-number distribu- 
tion for certain values of squeeze parameter. These oscillations are purely quantum effect 
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FIG. 2: Photon-number distribution P{n\) for the signal mode for oy = 3exp(i?), j = 1, 2, 3, £ 
1.5 and Ai = 0.3. The dashed curve is the Poisson photon-number distribution. 



without classical analogue and they can be interpreted as interference in phase space 



12] 



or in the framework of the generalized superposition of coherent fields and quantum noise 



351 ] . Indeed, these oscillations can be recognized for the linear mode where this mode can 
display squeezed light (see Fig. 3, for shown values of parameters). From this figure we see 
the macroscopic oscillations of the photon- number distribution for squeezed light. Further, 
comparison of Figs. 3a and 3b shows that the nonclassical oscillations of P{n%) are faster 
and more pronounced when t increases. Moreover, we noted that these oscillations appear 
only when the intensities of the input light are weak otherwise the single-peak structure is 
dominant. The broader oscillation range means that {(Ah^t)) 2 ) > (n 3 (t)) and thus there 
are always super-Poissonian statistics for the linear mode. It is worthwhile mentioning that 
the case of Fig. 3a corresponds to squeezed light with squeeze parameter r ~ 3.6. 

In fact, the behaviour of the photon- number distribution in this device is slightly different 
from that of the nonlinear asymmetric coupler with strong classical stimulating light in the 
second harmonic mode [23) where oscillatory behaviour as well as sub-Poissonian statistics 
for specific modes in both linear and nonlinear waveguide are exhibited. 
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FIG. 3: Photon-number distribution P{n^) for the linear mode when ay = 0.5exp(i^),j 
1, 2, 3; Ai = 1 and a) t = t s ; b) t = 5. 

V. PHASE DISTRIBUTION 



Nonlinear optical phenomena are sources of optical fields, the statistical properties of 
which are changed as a result of the nonlinear interaction. Quantum phase properties 
are among those statistical properties which undergo nonlinear changes. The interest in 
the phase properties has been motivated by experimental realization of optical homodyne 
tomography 36] allowing quantum phase from the measured field density matrix. 

We make use of the single-mode Q-function to investigate the phase distribution for the 
system under discussion by integrating this function over the radial variable 



P(Q,t) 



W(P,s = -l,t)\P\d\P\- 



It is worth mentioning that the authors of [37 



(23) 
have defined a Hermitian relative phase 



operator, which is in fact a polar decomposition of the Stokes operators for the two-mode 
field. In order to measure the single-mode phase distribution in this approach we have to 



prepare one of the modes in a very intense state o: 
information of the other mode can be estimated 



well defined phase and then the phase 



experimentally demonstrated 



3- 



Recently, this approach has been 



Actually, this approach gives different behaviour than 



that proposed from quasiprobability distribution functions for weak fields because of variety 
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of possible definitions of quantum phase 131 ] . However, they are in a good agreement for 
strong fields 40]. 



We proceed, equations (13) and (23) lead to 



7rC(t)v/A + (t)A_(t) 




cxp 



Ucw; 



2[A_(t) + l] 
+ erf(- 



2[A+(t) + 1] 



L4 + (i)sin 2 + yL(£)cos 2 0], 

{A.(*)[ai(t) + a\{t)} sin6 - M+(t)[ai(t) - a\{t)) cos0} 



(24) 



A_(t)i4+(i) 

and we have used the Gauss error function, which is defined by the well-known formula 

2 r 

erf (a) = —= I exp(-y )dy, 
V 71 " Jo 
taking into account that /3 = |/3| exp(iO) 



(25) 

i 

(26) 



It is clear that (24) is 27r-periodic function and the behaviour of this distribution in this 
device can be understood on the basis of the competition between the two-peak structure 
for the vacuum states and the single-peak structure for coherent states provided that the 
coherent amplitudes are real. Such behaviour is demonstrated in Figs. 4a,b for the signal 
mode and in Figs. 5a,b for the linear mode, where the evolution of the phase distribution 
(24) has been depicted against the phase and the time t. 

In figures a and b, oij = and 1, j = 1,2,3, respectively, for the shown value of coupling 
constant Ai. In Fig. 4a, one can see the phase evolution for the vacuum states in the signal 
mode showing two peak-structure with one central peak at = and two wings as — > ±7r. 
In fact, this behaviour has been seen for some of superposition states, e.g. for Yurke-Stoler 
states 41] and odd-binomial states [42J]. In Fig. 4b, one can observe that the initial peak for 



coherent state is attenuated in the course of time while the two peaks for the vacuum states 
are created and amplified. The well-known behaviour for the phase distribution of squeezed 
states may be established when we focus our attention on the behaviour of the linear mode, 
as expected (see Figs. 5a and b). 

From figure 5a one can see that the two-peak structure is typical for squeezed vacuum 
states, i.e. we have two-peak structure for = ±| 13] . Further, the heights of these 
peaks are j^\ a (t) a ^ an y time t (t > 0). However, Fig. 5b displays the well-known 
bifurcation shape for the phase distribution of squeezed field, i.e. the distribution curve 
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Phase distribution P(@, t) for the signal mode against and t for Ai = 0.6 and a) aj = 0; 
1, j = 1,2,3. 



undergoes a transition from single- to a double-peaked form with increasing time. Indeed 

n 

this figure is quite similar to that for squeezed coherent states [43], and the distributions 
differ in the behaviour of the initial peak, which is here amplified for a while before splitting 
into two peaks. This is connected with the significant influence of the development by 
nonlinear effects and power transfers between waveguides. It is worthwhile mentioning that 
a similar behaviour has been obtained for a contradirectional nonlinear asymmetric coupler 
with strong stimulated coherent field in the second harmonic waveguide [24| . 

Finally, we would like to conclude this section by discussing the origin of such behaviour 
of the phase distribution. This can be easily understood by analysing the function P(Q,t) 
when aj = 0, j = 1, 2, 3. In this case the formula (24) reduces to 



P(Q,t) 



^/A + {t)A_{t) 



(27) 



2vr A + (t) sin 2 9 + A_ (£) cos 2 9 ' 
This distribution function is double- or three-peak structure according to the relation be- 
tween A + (t) and A_(£). Let us restrict our discussion to the one of the modes which are 
propagating in the nonlinear waveguide, say, to the signal mode. For this mode, it can easily 
be shown that A+(t) > A-(t) and consequently the formula (27) exhibits three-peak struc- 
ture with peaks for 9 = 0, ±7r. The heights of these peaks at any time t ^ are -L< / A + \1 - 
For non-zero displacements aj, an additional factor of the form for coherent state, but with 
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FIG. 5: Phase distribution P(@,t) for the linear mode against and t for the same situation as 
in Fig. 4. 

a coherent amplitude ai(t) which is the expectation value of the signal mode operator in 
coherent state, appears in the distribution. This factor is responsible for a peak at = 
related to coherent component, which competes with the three-peak structure of vacuum 
to provide the previous behaviour of the distribution. Similar argument can be adopted to 
explain the behaviour of the linear mode, however, it should be borne in mind that in this 
case A + (t) < A-(t), as we have shown earlier for squeezing phenomenon, and hence the 
two-peak structure for the input vacuum case is dominant at 6 = ±|. 

VI. CONCLUSION 



In this article we have investigated quantum statistical properties of squeezed light gen- 
erated in a nonlinear asymmetric directional coupler composed of a linear waveguide and a 
nonlinear waveguide operating by nondegenerate parametric amplification. The main dif- 
ference between the present system and the well documented optical parametric oscillator 
(OPO) is the operating mechanism of these devices. More illustratively, OPO is using 
a quadratic nonlinear medium [9[ in a resonator and it is a practical device to generate 
squeezing; the system discussed here is richer by construction and by their properties be- 
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cause it includes two propagating media (i.e. waveguides), where one is linear and the other 
is nonlinear and they exchange energy linearly via evanescent waves. The system consid- 
ered here consists of a mixture of nonlinear processes including parametric generation and 
amplification (interaction coefficient Ai) and frequency conversion (interaction coefficients 
A2 and A3). This can lead to a single mode generation of nonclassical light. 

Using the Heisenberg approach we have examined quantum statistics of interacting 
modes, including the linear mode, propagating in the linear waveguide, in terms of the 
quasiprobability functions, photon-number distribution and phase distribution when the 
modes are initially in coherent states. From the behaviour of the Glauber P- and Wigner 
VU-functions in the linear mode it follows that the former is not well-behaved function when 
t > 0, while, the latter displays the well-known behaviour for squeezed light in which one of 
the quadratures is amplified and the other is attenuated. In the photon-number distribution, 
the large scale macroscopic oscillations related to squeezed light are established. The phase 
distribution displays the bifurcation typical for squeezed field. This behaviour demonstrates 
that changing the power of the linear interaction, it is possible to transfer nonclassical prop- 
erties, which are generated in the nonlinear waveguide, to the linear signal mode and to 
control them. Further, we have shown, in general, that the system is more suitable for 
generation of squeezed light rather than sub-Poissonian light from initial coherent light. 
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